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The magnetic state of Nitrogen-doped MgO, with N substituting O at concentrations between 1% 
and the concentrated limit, is calculated with density-functional methods. The N atoms are found 
to be magnetic with a moment of 1 per Nitrogen atom and to interact ferromagnetically via the 
double exchange mechanism. The long-range magnetic order is established above a finite concentra- 
tion of about 1.5% when the percolation threshold is reached. The Curie temperature Tc increases 
linearly with the concentration, and is found to be about 30 K for f0% concentration. Besides the 
substitution of single Nitrogen atoms, also interstitial Nitrogen atoms, clusters of Nitrogen atoms 
and their structural relaxation on the magnetism are discussed. Possible scenarios of engineering 
a higher Curie temperature are analyzed, with the conclusion that an increase of Tc is dicult to 
achieve, requiring a particular attention to the choice of chemistry. 

PACS numbers: 75.50.Pp, 75.50.Hx, 75.30.Et 



I. INTRODUCTION 

In the research field of diluted magnetic semiconduc- 
tors (DMS), a new direction is being investigated in the 
last five years, namely the engineering of ferromagnetic 
state formation by sp impurity doping. Compared to 
the more traditional DMS with transition-metal impuri- 
ties, the novel sp-magnctism, or d°-magnctism, is rather 
unexplored. The increasing interest is due to the per- 
haps unexpected finding that p-bands can spontaneously 
polarize giving a ferromagnetic state (although the pos- 
sibility of magnetic sp-defects is long known), but also 
due to the hope of tuning the properties of these states 
so that high Curie temperatures are achieved even at low 
concentrations . 

Theoretical considerations have revealed mainly two 
routes for the formation of sp-ferromagnetic states. In 
the first, the semiconductor or insulator cation is substi- 
tuted by an atom of smaller valency, thus depriving the 
p-typc valence band from electrons. This hole-doping can 
shift the Fermi level into the valence band deep enough 
that the Stoner criterion is fulfilled, and spontaneous spin 
polarization appears (see Fig. [TJi). This is, for example, 
the mechanism encountered in alkali-atom doped T1O2 
and Zr02, predicted to be ferromagnetic by ab-initio 
calculations i The hole-doping can also be achieved by 
cation vacancies, instead of substitution, as was e.g. pro- 
posed for the cases of Hf02, CaO, and ZrC^^^ 

The second scenario is that the anion is substituted 
by one of smaller valency, introducing shallow, spin- 
polarized gap states. As the impurity concentration 
increases, these states form impurity bands, which re- 
main spin-polarized if the Stoner criterion is fulfilled 
(see Fig. [Tp). The magnetic moment is expected to be 
strongest for 2p impurities, i.e., in the case in Carbon- or 
Nitrogen-doped oxides^&^&^ i 10 ' 1 ^ 12 ' 13 This is because 
the 2p states have no nodes and are rather localized, 
leading to a significant Hund's-type exchange interaction 
on-site. In the concentrated limit (full substitution of the 
anion, see e.g. Ref. [13), both routes converge to the same 
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FIG. f : (color online) Typical mechanisms for sp ferromag- 
netism: (a) cation substitution drives the Fermi level _Ep into 
the valence band; (b) anion substitution creates an impurity 
band in the gap. 



mechanism. 

While the appearance of p-type magnetic states in ox- 
ides has been studied in the past, preliminary predic- 
tions of the Curie temperature have been so-far based 
only on mean-field theory^&ii which has since been 
shown to give qualitatively wrong results, overestimat- 
ing Tc, because it ignores the fundamental phenomenon 
of magnetic percolation in diluted disordered magnetic 



systems 



15.16.17 



This means that, at low concentrations, 



the average inter-impurity distance is large, so that the 
short-range exchange interactions cannot produce a high 
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Curie temperature although they are strong. 

The scope of this paper is to investigate the appear- 
ance and stability of the ferromagnetic state, including 
calculations for the Curie temperature, in MgOi-^N^ 
compounds for x < 15%, in view of recent experimen- 
tal activities in the particular system.— Our focus is on 
the solution provided by local density-functional theory 
(DFT), presented in Sec. UH and on the calculation of 
the Curie temperature (Sec. IIII|) with the exchange con- 
stants harvested within the adiabatic approximation. As 
it turns out, for any reasonable concentration the Curie 
temperatures are significantly lower than room tempera- 
ture, therefore, we also explore the possibility of increas- 
ing the Curie temperature by heavy doping in Sec. IIVI 
Further on, we discuss the physics that lies beyond our 
approximations in Sec. [V] There, among our consider- 
ations, a comparison is made to previously calculated 
results for N-doped MgO that focus on the effects of 
electron correlation, and we also explore the constraints 
imposed by our specific structural model, including a dis- 
cussion on the solubility of N in MgO. Finally, we give 
an outlook in Scc lVIl The methods of calculation arc 
shortly described in the Appendix. 



II. GROUND-STATE ELECTRONIC AND 
MAGNETIC STRUCTURE 

MgO crystallizes in the rock-salt structure with a lat- 
tice parameter of 4.21 A and exhibits a wide band gap 
of 7.8 cV. In our density-functional calculations the gap 
is found to be 4.8 eV, due to the well-known under- 
estimation of insulator band gaps within local density- 
functional theory [local density approximation (LDA) or 
generalized gradient approximation (GGA)]. When sub- 
stituting O, N induces p states in the MgO gap, approx- 
imately 0.5 eV above the valence band edge. As N has 
one electron less than O, the gap states host one unpaired 
hole, showing a magnetic moment of 1 /ib per N atom. 

At finite concentrations, interaction among the impu- 
rity states forms a partially filled impurity band. The 
spin polarization remains present even at high concen- 
trations, showing that magnetism is not only a conse- 
quence of the missing electron, but also of the relatively 
strong Hund's exchange due to the localization of the 
2p states. The spin-polarized density of states shows a 
half-metallic behavior, as we found by calculations with 
the Korringa-Kohn-Rostoker Green function method^! 
(KKR) within the Coherent Potential Approximation to 
disorder (CPA). (Details on the calculation methods are 
given in the Appendix.) Through a calculation of struc- 
tural relaxation around a N impurity, performed with 
the VASP— projector augmented wave code at 5% con- 
centration in supercell geometry, we could rule out seri- 
ous Jahn- Teller distortions around the impurity: the first 
neighbors relax outward by less that 2%, including only a 
weak symmetry breaking along one of the three crystallo- 
graphic axes. This leads to a small lifting of degeneracy 
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FIG. 2: (color online) (a) Atom-resolved local densities of 
states (LDOS) of MgO0.95N0.05 calculated within the KKR- 
CPA. The LDOS of the Mg atom is presented by a dashed 
line, (b) Same for MgO0.45S0.5N0.05- 



of the impurity state, insignificant compared to the im- 
purity band width. The small atomic displacements also 
justify the use of CPA. 

We turn now to the discussion of the electronic and 
magnetic structure. Fig. [2ja) shows the atom-resolved 
local densities of states (DOS) at a N concentration of 
x=5%. Evidently N has a spin-split DOS and exhibits a 
local magnetic moment. The N majority-spin (spin-up) 
impurity band is fully occupied and the minority-spin 
(spin-down) impurity band is occupied by 2/3. As a con- 
sequence, 1/3 the total moment is 1 /ib per N atom. 
At lower concentrations the DOS is very similar, but 
with smaller impurity band width w (w depends on the 
concentration as w ~ y/x, since x represents an aver- 
age number of impurity neighbors, and it is known from 
the tight-binding approximation that w increases as the 
square root of the number of neighbors). The exchange 
splitting is of the order of 0.5 eV, giving an exchange in- 
tegral of / = 0.5 cV/yUB, which is approximately half the 
one in 3d transition metals. 

Already from examining the ferromagnetic density 
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FIG. 3: Spin-polarization energy per Nitrogen atom as a func- 
tion of concentration x. As the concentration increases, the 
hybridization increases, the impurity band becomes more itin- 
erant, and the spin-polarization energy per atom drops. The 
line is a guide to the eye. 



of states it is expected that the ferromagnetic state 
will be stable, since E-p lies within the impurity band, 
favoring double-exchang e 20 ' 21 (i.e., upon ferromagnetic 
alignment, hybridization leads to a broadening of the 
partially-filled impurity-band, resulting in energy gain). 
We confirmed this by calculations (not shown here in de- 
tail) of the total energy of the ferromagnetic state, -Ef er ro, 
versus the disordered-local- moment state energy, £dlm; 
the latter is represented in the CPA by an alloy of the 
type MgOi_ x Nj. 5a .N^ Sa ., where N T and are impuri- 
ties with magnetic moment pointing "up" and "down" 



As is typical for double-exchange DMS systems,— we 
found that the energy gain of the ferromagnetic state, 
-Edlm — Efe ITO , scales with the square root of the concen- 
tration. Further confirmation on the ferromagnetic na- 
ture of the ground state comes from calculations of the 
exchange constants, to be presented in the next Section. 

The spin polarization energy (i.e., the energy gain of 
the system due to the moments' formation) drops with 
concentration due to the increase of hybridization, rang- 
ing between 550 and 100 meV for 2% < x < 15%, as 
can be seen in Fig. [3J However, a high spin-polarization 
energy is a necessary but not sufficient condition for the 
stability of the ferromagnetic state and for a high Curie 
temperature. The extension of the impurity gap-state 
is also important, as it dictates the "communication" 
between localized moments and determines the onset of 
magnetic percolation i 15 ' 16 We therefore show in Fig. [D(a) 
the extension of the Nitrogen-induced spin-polarized hole 
(calculations here were done within the KKR impurity- 
in-host approach, assuming that the three N-induced gap 
states are degenerate and occupied by 2/3 each). We 
find that approximately half of the hole is localized at 
the impurity site, while most of the remaining charge is 
distributed at the 12 nearest Oxygen neighbors. From 
Fig. 0Jb) it can also be seen that the hole falls off ex- 
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FIG. 4: (color online) Spatial distribution of the Nitrogen- 
induced spin-polarized hole for the case of a single Nitro- 
gen impurity substituting O in MgO. Only the fraction of 
the hole at the N and O atoms is presented (the fraction at 
the Mg atoms is negligible). The distribution of the mag- 
netic moment (1 /xb) is practically the same, (a) Full (green) 
curve: Integrated hole charge as a function of distance from 
the impurity. Dashed (red) curve: Hole charge per shell of 
O atoms that are equidistant from N and equivalent by the 
point-group symmetry. The impurity accommodates more 
than half of the hole charge, (b) Hole charge per atom as 
a function of distance from the impurity. Only one of the 
symmetry-equivalent atoms is considered for every shell. The 
points labelled (110) represent atoms in the (110) direction, 
in which the extent of the impurity wavefunction is largest, 
ao is the lattice parameter. 



ponentially with distance, showing, however, a higher 
value in the (110) directions. These are the directions 
of the first Oxygen neighbors in the fee lattice (which is 
the sublatice formed by Oxygen atoms in MgO), and in 
these directions also the pair exchange constants fall off 
slower (see Sec. IIII|) . The hole occupation at Mg sites is 
negligible. 
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FIG. 5: (color online) Pair exchange interactions J nn i as a 
function of distance R between nitrogen atoms for x = 1%, 
3%, and 5%. The peaks corresponding to neighbors in the 
(110) directions are indicated, together with the percolation 
threshold for interactions up to the particular distance, ao is 
the lattice parameter. The lines are guides to the eye. 



III. EXCHANGE INTERACTIONS AND CURIE 
TEMPERATURE 

We turn now to the discussion of the Curie temper- 
ature. We describe the fluctuations of the magnetiza- 
tion on the basis of the classical Heisenberg Hamilto- 
nian, H = - J2 nn ' Jnn> e„ • e n >. Here, J nn > are the pair 
exchange parameters between N impurities at sites n, n' 
and e„, ev are unit vectors pointing in the direction of 
the local moments. The Heisenberg model includes the 
transversal degrees of freedom of the fluctuating mag- 
netic moments. As we see below, this energy scale is 
much smaller than the spin-polarization energy, and thus 
the longitudinal fluctuations of the magnetic moments 
are safely ignored as regards the Curie point. 

The J nn i are fitted to LDA results so that at the end 
the magnetic excitation energies of the Heisenberg Hamil- 
tonian correspond to the ones of MgOi-^Na;. For the fit- 
ting we employ the method of infinitesimal rotations jS 2 - 
The results are shown in Fig. \5\ Clearly, all interactions 
are ferromagnetic. We note that there is a marked con- 



centration dependence, with stronger exchange constants 
at low concentrations, as is well-known in DMS systems 
with double-exchange ferromagnetic interactions (see, for 
example, Ref . |23T) . 

Note, however, that there is a "hidden" antiferromag- 
netic superexchange interaction, reducing the value of 
Jnri ■ This arises from the proximity of the impurity 
bands of both spins to E-p: in an antifcrromagnctic 
alignment of two impurities, energy is gained due to hy- 
bridization of majority-spin states of each impurity with 
minority-spin states of the other, accompanied by down- 
ward shifting of the occupied majority-spin bands and 
upward shifting of the minority-spin bands^ This su- 
perexchange counter-acts the double exchange to some 
extent. To reveal this, we performed a non-self-consistcnt 
calculation of J nn i (just one iteration) after shifting the 
majority-spin potential of N to lower energies by 2.7 cV. 
In this way the majority-spin N band is driven away from 
the Fermi energy and superexchange is strongly reduced. 
The result (not shown here) was striking: the J nn i be- 
came stronger by approximately 50% even at long dis- 
tances. 

It is evident that the nearest-neighbor exchange con- 
stants are large (Fig. O, though insignificant, in the di- 
lute case, for the stability of the ferromagnetic state; 
this is governed by longer-range exchange, due to the 
requirement for magnetic percolation. 15 ' 16 However, the 
J nn i fall off exponentially with distance, as is expected 
by the fact that the Fermi level falls in the gap at least 
for the one spin directional (half-metallic or insulating 
systems). Notably, the slowest decay of J nn i is observed 
along (110), i.e., in the directions where the charge of the 
impurity gap-state falls off with the slowest rate (see also 
Fig. [4]). Note that the interactions in the (110) directions 
are also dominating in zinc-blende or diamond-structure 
DMS with transition-element impurities] 25 ' 26 

In Fig. Owe also give the percolation thresholds which 
we calculated for interaction distances corresponding to 
"peaks" of J nn i in the (110) directions. We see that the 
nearest-neighbor coupling starts playing a role at x = 
20%, while the next peak, corresponding to R = s/2ao, 
becomes important only for concentrations above x = 
4.9%; at this concentration, J(v / 2a ) ~ 5 meV. At con- 
centrations of the order of 1.5%, where we see the onset 
of the Curie temperature in Fig. [5] (discussed in more de- 
tail below), the third neighbor in the (110) direction be- 
comes important, with a low value of J(^ ao) ~ 1 meV. 
We consider interactions beyond this distance as negli- 
gibly small, therefore we recognize the concentration of 
1.5% as a magnetic percolation threshold for MgOi-zN^; 
this can be seen also from the concentration-dependent 
Curie temperature, as we discuss below. (Of course J(R) 
is never exactly zero, therefore the choice of x = 1.5% 
as percolation threshold is somewhat arbitrary, but, we 
think, reasonable.) 

A strong exponential decrease of J nn i is an indica- 
tion of a low Curie temperature in diluted magnetic sys- 
tems. We have calculated Tq within the random-phase 
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FIG. 6: Curie temperatures for N concentrations between 1% 
and 15%, calculated within the RPA. Inset: Same in the con- 
centration range below 1-90%. The lines are guides to the 
eye. 



approximation for disordered systems, following the pre- 
scription of Hilbert and Noltingii (see the Appendix for 
calculational details). In doing this, we assumed the 
limit of classical spins (S — ► oo, with appropriate re- 
normalization of the exchange constants), for reasons 
that we discuss in Sec. [V] In special cases we also per- 
formed Monte Carlo calculations, as described in the Ap- 
pendix, which gave typically 50% higher Tq. The results 
on the Curie temperature, as a function of concentration, 
are shown in Fig. [BJ Although in this work we are inter- 
ested in concentrations up to the order of x = 10-15%, 
in the calculation we included also high concentrations, 
in order to see the trend when approaching the concen- 
trated limit. Starting from a percolation threshold at 
approximately x = 1.5%, Tq rises linearly with concen- 
tration. Although, as we saw, the exchange constants 
at any particular distance weaken with increasing con- 
centration, the important effect here comes from the de- 
crease of the average inter-impurity distance at higher x, 
allowing for the shorter-ranged, stronger interactions to 
play a role in the ferromagnetic ordering. 

At concentrations that are usually considered in DMS 
systems, i.e., not higher than 10-15%, the Curie temper- 
ature turns out to be rather low: we find Tc = 35 K 
at 10%, and even at x = 20%, T c is still below 100 K. 
This is a clear drawback for application purposes at room 
temperature, and in the following Section we consider 
possible ways to overcome this difficulty. 

IV. CONSIDERATIONS ON HOW TO 
INCREASE T c 

Since the Curie temperature of diluted magnetic sys- 
tems depends on the long-range exchange constants, we 
have examined several possibilities of increasing J„„< at 
distances larger than nearest neighbors. Considering 



that J nn > falls off exponentially with distance R nn ' = 
\R n — Rn>\ with a characteristic decay parameter k, 
Jnn' ~ exp(— nRnn'), the main scope is to reduce k. We 
tried to do this in several ways, as we discuss here, how- 
ever, none produced a significant improvement. 

-Reducing the band gap. Perhaps the most obvious 
way to reduce k is to engineer a smaller gap by alloying. 
For instance, Mgi-^Zn^O shows a smaller gap, while re- 
taining the rock-salt structure at not too high Zn con- 
centrations; numerous other alloying combinations could 
have a similar effect. We pursued this idea by acting with 
an attractive constant potential on the Mg site, thus ar- 
tificially reducing the MgO gap size. We then derived 
k(E) by calculating the complex band structure for en- 
ergies in the gap region, i.e., looking for solutions of the 
Schrodingcr equation that fall off as ip ~ e~ KT (this pro- 
cedure is well-known in the calculation of surface and 
interface states). For E deep in the gap there is a signif- 
icant reduction of k(E) with the gap size. However, for 
E at the actual position of the N levels (relatively close 
to the valence band edge E v ), k(E) is mainly determined 
by the effective mass m* of the light-hole valence band 
asKRi y / 2m*(E — E v )/h, due to the analytical behavior 
of the complex band structure. (Note that, among the 
complex bands, we seek the one with smallest value of 
K, which corresponds to the band derived from the light 
holes As the effective mass does not change apprecia- 
bly, we find that a reduction of the gap does not affect 
n(E-p) significantly, thus it cannot change the long-range 
exchange coupling. 

-Changing the effective mass by compression. After 
the previous considerations, the next obvious idea is to 
reduce the effective mass by compression, to be achieved, 
e.g., by epitaxial stress or strain. Calculations for a 5% 
reduction of the lattice constant show that only the ef- 
fective mass of the heavy holes is appreciably affected, 
while the light hole bands are not much altered. How- 
ever, what is important for a small k is the light-hole 
behavior; thus k(Ef) docs not change practically, and 
there is no increase in J„ n >. Actually for large distances 
we do obtain small changes in J rm / after compression, but 
towards antiferromagnetic behavior. This has probably 
two origins^ at smaller lattice constants, the band width 
increases and the exchange splitting becomes smaller, so 
that (i) the spin-down impurity band merges more with 
the valence band setting on a Ruderman-Kittel-Kasuya- 
Yosida (RKKY) behavior, and (ii) the antiferromagnetic 
superexchange, which is sensitive to the exchange split- 
ting, increases. 

-Alloying with Sulfur. Next we considered shifting the 
N impurity states closer to the valence band by proper 
alloying of MgO with a third component. In this re- 
spect, MgOi_ y Sj, seemed promising, because MgS also 
crystallizes in the rock-salt structure, with a larger lat- 
tice parameter (5.2 A) and a smaller, indirect band gap 
(experimental: 4.5 eV ; we find 2.25 cV within the LDA). 
Calculations in MgOi-^Sy (which were carried out by 
changing the lattice constant according to Vegard's law) 



6 



showed that the O valence states end up higher than the 
S valence states (see Fig. d^b)). This result is counter- 
intuitive if one takes into account only the electronegativ- 
ity of the elements. We interpret it by observing that the 
stronger localization of the O and N 2p states is responsi- 
ble for a stronger on-site Coulomb repulsion, so that filled 
2p levels end up higher in energy than the more extended 
S 3p levels. As a consequence, the valence band edge in 
MgOi-y-zSyNj; is again Oxygen-dominated, but the po- 
sition of the N impurity band is closer to the Oxygen va- 
lence band edge, and there is a stronger hybridization of 
the N with the O states compared to MgOi-^N^, as can 
be seen in Fig.[2la,b). Because of this, the pair exchange 
constants are more extended in MgOi—y-^SyN^. Calcu- 
lations of MgO0.45S0.5N0. 05 showed an enhancement of 
J(R) at large distances by up to 10-20% compared to 
MgO0.95N0.05- Thus the Curie temperature increases by 
a factor of the same order, remaining small. 



-Doping with holes. In principle, the double-exchange 
mechanism for fcrromagnctism is most effective when the 
impurity band of one spin is exactly half-filled, while for 
the opposite spin is completely full (or empty). Then 
the energy gain in the ferromagnetic state by band- 
broadening is maximized. In MgOi-^N^ this is not the 
case, as the N spin-down band is filled by 2/3. Therefore, 
doping with holes could help, as Ep will be shifted down- 
ward to the middle of the impurity band. Note that we 
are not trying to achieve the analogue of Zener's hole- 
mediated p-d exchange, but rather an improvement of 
the double-exchange mechanism. (Similar engineering 
has been proposed in transition- metal-doped DMS.— ) 
We investigate this effect in a Mg-poor compound i.e., 
Mgi-^Oi-^Na;. Each missing Mg atom adds two holes; 
thus the "ideal" Mg-vacancy concentration, resulting in 
exactly half-filling of the N spin-down band, is y = x/4. 
In our calculations we used x = 10% and y = 2.5%. 
As it turns out, the exchange interactions become more 
extended in space, however, the reason is that Ep now 
slightly enters the O valence band (since the valence band 
top is anyhow merged with the lower part of the im- 
purity band), and a Fermi surface is now available for 
the impurity interaction. As a result, a RKKY mecha- 
nism emerges, with an outcome of also antiferromagnctic 
interactions at large distances. The latter turn out to 
be comparatively strong, once more because of an addi- 
tional increase of superexchange, as Ep shifts closer to 
the spin- up band. Thus there is no increase of Tc- Tests 
with lower hole-doping concentration, before the RKKY 
exchange sets in, result in only a marginal increase of 
the long-range interactions, insufficient for a significant 
increase of Tq, 



V. DISCUSSION ON APPROXIMATIONS 

A. Approximations in the calculation of the 
electronic structure and Curie temperature 

We now consider which approximations used in our cal- 
culations could most seriously affect the electronic struc- 
ture, exchange constants and Curie temperature. Per- 
haps the most critical is the local density approximation. 
Corrections that could qualitatively change the electronic 
structure include self-interaction correction, inclusion of 
strong on-site Coulomb repulsion, or GW-type of self- 
energy. 

In Ref. H3, Droghetti and Sanvito applied the self- 
interaction correction (SIC) to MgO doped with N 
(among other compounds). They find that an applica- 
tion of SIC to both N and O atoms increases the gap, 
as is expected by driving the occupied levels lower by 
the SIC. However, their main finding connected to our 
present discussion is the occurrence of a strong splitting 
of the order of 3 eV between the occupied and unoccupied 
p states of N. The spin-down impurity band is, after SIC, 
manifestly insulating. A similar effect was found within 
the LDA+/7 method, applied by Pardo and Pickett to 
a number of oxides (including MgO) with N substitution 
of Oxygen. In the latter work, a Coulomb repulsion of 
U = 5.5 eV was used on the O and N atoms. From the 
point of view of ferromagnetism, an insulating impurity 
band can have severe consequences, in particular disfa- 
voring the double-exchange mechanism (impurity-band 
broadening docs not lead to energy gain any more) ; ferro- 
magnetic superexchange could be present, but it is much 
weaker. In these cases, the system could even show a 
spin-glass ground state instead of a ferromagnetic one. 

These works demonstrate possible consequences of 
strong Coulomb interactions, showing that local density- 
functional theory could be insufficient. We believe that, 
at low concentrations, these results undoubtedly show 
the correct physics, but at higher concentrations (a few 
percent) they possibly overestimate the Coulomb inter- 
action, and are therefore still inconclusive as to the exact 
nature of the ground state. Our arguments are as follows. 

First, considering the low-concentration limit, the N 
impurity states are not localized at the N atomic cell, 
but rather extended also over the 12 nearest Oxygen 
neighbors, as our calculation shows. This speaks for 
a relatively mild self-interaction error, compared to d- 
or /-systems, where the SIC is usually applied to. The 
LDA error must become less and less serious as the con- 
centration increases and an itinerant impurity band is 
formed. Furthermore, application of Coulomb correction 
terms on the Oxygen valence band is probably much less 
necessary, as the associated Bloch states arc itinerant 
with a significant band width (in the limiting case of 
completely non-localized states of a homogeneous elec- 
tron gas, the self-interaction is fully corrected by the 
LDA exchange-correlation energy). We conclude that 
Coulomb-interaction corrections should be of different 
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scale at the N impurity-band states compared to the O 
valence-band states, but even for N they should not be 
too severe. 

Second, as the concentration x increases, the impurity 
band width w ~ \fx increases rather fast, exceeding 1 eV 
already at x = 5%. This band width can well be of the 
order of magnitude of the local Coulomb interaction U (in 
the LDA+C/ calculations of Ref . [Til, a value of U = 5.5 eV 
was used, which, however, acts on the atomic site; if the 
full extent of the impurity wavcfunction is accounted for, 
a significantly lower value of U would be needed for the 
same effect). In this case a Mott transition is possible, 
from an insulating state for low concentrations (where 
w < U) to a metallic state at high concentrations (where 
w > U). In such a scenario, ferromagnetism would be 
assisted by the increase of concentration in two ways: 
occurrence of a metallic state similar to the LDA result, 
and magnetic percolation. 

Concerning GW-type of corrections, these are known 
to correct the LDA underestimation of band gap in band 
insulators. It should be then expected that the larger gap 
would lead to a stronger decay of the exchange constants 
with distance, reducing the calculated To However, this 
reduction should not be too serious, for the same rea- 
son that Tc cannot be much increased by engineering a 
smaller gap (see Sec. IIV[) . 

Another approximation that was made here was the 
assumption of a classical, rather than quantum, Heisen- 
berg model to describe the magnetic excitations. Con- 
sidering that we are faced with an S = | system, it can 
be argued that a classical approximation is unrealistic. 
Within the RPA, a change from a classical to a quan- 
tum Hcisenberg model leads to an increase of Tc by a 
factor of S(S + l)/S 2 = 3, in the S = \ case. Ac- 
cepting this, room-temperature ferromagnetism could be 
achieved at x = 20%, as can be seen by scaling up the 
results of Fig. [6l However, the particular way of calcu- 
lation of the exchange coupling constants^ J nn i tacitly 
assumes a classical model. The constants are calculated 
within constrained density-functional theory, in principle 
by "freezing" the system in a static non-collinear config- 
uration and calculating the total energy; this should be 
viewed as a parametrization of the low-energy excitations 
of the spin density, rather than a derivation of a quantum 
spin Hamiltonian. Comparison to experiments in previ- 
ous works also advocates for this point of view. E.g., 
Sasioglu et alM- have calculated the Curie temperature of 
Heusler alloys by a similar recipe. While the assumption 
of a classical Heisenberg model lead them to reasonable 
agreement of Tc with experiment, the assumption of a 
quantum model resulted in a clear overestimation of Tc . 
For lack of a better theory that predicts J nn i for a quan- 
tum Heisenberg model, we are obliged to work within a 
classical model. 



B. Approximations imposed by our structural 
model 

So far we assumed a uniform, on the average, distri- 
bution of the N atoms in the MgO matrix. We there- 
fore start this Subsection by commenting on the sol- 
ubility of N in MgO. The solubility limit of substitu- 
tional N is expected to be small, and clustering of the 
N impurities is expected to be favored. We verified 
this by calculating, within the KKR-CPA, the mixing 
energy as a function of the concentration, E m { x (x) = 
TpVIgO^Rr] - (1 - x)T[MgO] - xE[MgN] with MgN 
in the rock-salt structure. E m { x (x) was positive for all 
calculated concentrations (1% < x < 90%), showing 
that phase separation into MgO and MgN is energeti- 
cally favored. This can be understood also from the lo- 
cal density of states, shown in Fig. [Ha), where it is seen 
that the impurity band is bisected by the Fermi level 
Tp. In case of impurity clustering the bandwidth will 
increase and energy will be gained because of the low- 
ering of the occupied levels. Interestingly, the origin of 
ferromagnetic double exchange (bisection of the impurity 
band by Ep) provides also a mechanism for (usually un- 
wanted) phase separation. Of course, phase separation 
other than the MgO-MgN type could be also possible; 
what this calculation shows is the thermodynamic insta- 
bility of MgOi-^Nj. with respect to at least one type of 
phase separation, even if entropic effects could milden 
the separation. The preference toward clustering was 
also seen in calculations of two substitutional N impu- 
rities in a MgO supercell matrix: the most stable state 
was found when the two N atoms were first neighbors in 
the (110) direction. We conclude that MgOi_ x N x can 
only be grown under out-of-equilibrium conditions, as is 
the case with many transition-metal-doped DMS. 

Two neighboring substitutional N atoms still provide 
a metallic state. However, three neighboring N atoms 
forming an equilateral triangle in a {111} plane cause 
a transition to a semiconducting state (with a split 
minority-spin band), if structural optimization of the 
three atoms and their surroundings is taken into account. 

Coming now to interstitial N impurities, our calcula- 
tions (performed with VASP) showed that strong atom 
displacements are involved, thus we cannot use CPA to 
describe the random alloy with interstitials. As the phase 
space involved is extremely rich, and probably depends 
strongly on growth and annealing conditions, we defer 
a deeper analysis of the particular question to a future 
work, commenting on a few interesting findings. 

First, in the case of a single N interstitial, the symmet- 
ric, tetrahcdral position constitutes only a local minimum 
of the total energy. Further significant reduction of the 
total energy is found in a "dumbbell" configuration rem- 
iniscent of a NO molecule, where the N impurity binds 
itself to an O atom, with both atoms in the proximity of 
the ideal O lattice position. This entity is found to be 
half-metallic and magnetic, with a moment of 1 [1b- 

Second, when two Nitrogen atoms are placed around 
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an Oxygen atom as interstitials, structural relaxation 
leads to a "zig-zag" O-N-O-N configuration. Here, the 
two N atoms are almost on the lattice sites, being first 
neighbors along (110), while the O atoms are shifted out 
of their ideal positions, hovering above the (110) atomic 
layer. The electronic structure changes to non magnetic 
and insulating. 

Third, the following configuration is of particular in- 
terest: two neighboring N atoms together with an O in- 
terstitial between them, plus a nearby O vacancy. This 
is so to say a configuration where an O atom has been 
bound during growth by a N pair, missing the lattice 
position nearby. This configuration shows a local en- 
ergy minimum and is non-magnetic. However, if the O 
atom returns to the vacant lattice position, the state be- 
comes magnetic with a total energy gain of the order of 
4 cV, while the lattice parameter is reduced by 1% (in 
the particular calculation we took a Mg32C>3oN2 super- 
cell, corresponding to 6.25% N). This is a possible expla- 
nation of the appearance of magnetism after annealing 
in experiment^ 



VI. CONCLUSIONS AND OUTLOOK 

The chemical compound MgOi-^N^ in which Oxygen 
is substituted by Nitrogen atoms is found to be ferro- 
magnetic for all concentrations which seem realistically 
achievable (x < 15%) and above the concentration of the 
percolation threshold (about 1.5%). The N atom forms 
a partially occupied gap state at the vicinity of the va- 
lence band edge. This leads to concentration-dependent 
exchange parameters, exponentially decreasing with dis- 
tance, which are ferromagnetic due to double exchange as 
dominating mechanism. The reduction of the average N- 
N distance and the reduction of the exchange interaction 
strength with increasing concentration compensate par- 
tially, so that the Curie temperature increases linearly 
with concentration. The structural relaxation around 
the N impurity are found to be small. In the case of N 
clusters and interstitial N we found magnetic as well as 
non-magnetic complexes, but never with higher moments 
than 1 ^B/atom. 

We think that the Curie temperature of such sp- 
compounds can be increased if the exchange interactions 
become more long-ranged, e.g. by moving the gap states 
closer to the band edge, and if the N-N exchange inter- 
actions are enhanced by suppressing the competing anti- 
ferromagnetic superexchange via an increase of the local 
exchange splitting. Our attempts in this direction, using 
compounds such as MgOi-^-ySyNj, or Mgi-yOi-^N^, 
have not been met with success. Nevertheless, it is worth- 
while to investigate more host/impurity combinations 
theoretically, in order to guide experimental efforts. 

From the point of view of applications, a low Tq is 
a limiting factor. However, from the point of view of 
physics, sp magnetism in oxides presents intriguing open 
problems, such as the nature of the ground state, its con- 



centration dependence, and the role of dynamic electron 
correlations. Moreover, although homogeneous Nitro- 
gen distributions at high concentrations are difficult to 
achieve, it is worthwhile to examine the physics and tech- 
nological relevance of inhomogeneous samples, emerging, 
e.g., by spinodal decomposition^ or delta-doping^ as 
has been proposed for transition- metal doped DMS. Such 
compounds could show much more robust magnetism lo- 
cally, leading to new functionalities. 



APPENDIX A: METHODS OF CALCULATION 

For the calculations of the alloy ground-state electronic 
structure (except structural relaxations) we used the 
full-potential Korringa-Kohn-Rostokcr Green-function 
method 3 -! (KKR) with exact calculation of the atomic cell 
shapes,— using an angular momentum cutoff of Z max = 3 
and a maximal fc-point mesh of 64000 points in the full 
Brillouin zone; the complex-energy contour integration 
included the Oxygen and Nitrogen 2s states. Disorder 
was treated within the Coherent Potential Approxima- 
tion (KKR-CPA). Rclativistic effects were neglected. In 
the KKR-CPA calculations we used the MgO experimen- 
tal lattice parameter for all concentrations, as the small 
change should have no significant effect in the calculated 
trends. 

For the calculation of the hole localization, the 
impurity-in-host version of the KKR method was em- 
ployed non-self-consistently, with potentials taken from 
self-consistent CPA. In the KKR-CPA calculations, the 
LDA for the exchange-correlation energy was used with 
the parametrization of Vosko ct ali^i For the struc- 
tural relaxations we used the Vienna ab-initio simulation 
package^ (VASP) within the GGA& 

The exchange constants J nn > were calculated within 
the approximation of infinitesimal rotations^ again 
within the KKR-CPA. As has been found by Sato et 
al.,— the method of infinitesimal rotations within the 
CPA can overestimate the values of J nn > between the 
nearest neighbors at low concentrations (when they any- 
how are irrelevant for Tq), but is rather good for larger 
distances. 

The Curie temperatures were calculated within the 
random-phase approximation (RPA) for disordered sys- 
tems, as proposed by Hilbert and Noltingjil with the 
distance-dependent interactions </„„< taken as described 
above. At each concentration, an environmental average 
was evaluated by taking 100 random configurations of ap- 
prox. 8500 nitrogen atoms each, statistically distributed 
in a simulation supcrcell (sized between 21 3 unit cells for 
90% concentration and 94 3 unit cells for 1% concentra- 
tion). To each configuration the RPA yielded a Tq, and 
an average of all values of Tq was determined at the end. 
In special cases (x = 5%, 10%, 15% and 10%) the results 
were cross-checked with Monte Carlo calculations, which 
yielded a Tc of approximately 50% higher, except for 5% 
where the results of the two methods agreed. Because 
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the Monte Carlo method is numerically more expensive, 
in these tests we treated systems of 1000-2500 magnetic 
ions, using up to 20 configurations for averaging; within 
Monte-Carlo, Tq was calculated using the 4th order cu- 
mulant method^ 
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